Fluctuations and isentropes near the chiral critical endpoint 
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Isentropic trajectories crossing the chiral phase transition near the critical endpoint (CEP) are 
studied for two light quark flavors. The calculations are performed within an effective chiral model 
with quark-meson interactions, belonging to the same universality class as QCD. We confront mean- 
field thermodynamics with the functional renormalization group approach, where fluctuations are 
properly taken into account. We establish a connection between modifications of the isentropic 
trajectories found in mean-field calculations at the crossover transition near the CEP and the order 
of the phase transition in the chiral limit. Furthermore, the isentropes obtained with the renormal- 
ization group are completely smooth at the crossover transition and do not in any way reflect the 
proximity of the CEP. In particular, our results do not show the recently conjectured focusing of 
isentropes from the crossover region towards the critical endpoint. 

PACS numbers: 



1. INTRODUCTION 

The determination of the phase structure and the crit- 
ical behavior of strongly interacting matter is one of the 
major goals addressed theoretically in studies of QCD at 
finite temperature T and quark chemical potential \x and 
experimentally in ultrarelativistic nuclear collisions. One 
of the fundamental predictions of QCD, is the existence 
of a boundary line in the (T, /z)-plane that separates the 
confined, chirally broken hadronic phase from the decon- 
fined chirally symmetric quark-gluon plasma (QGP) 

The existence of such a phase boundary was re- 
cently established by first-principle calculations in Lat- 
tice Gauge Theory (LGT) 0. For vanishing fi and for 
two massless flavors, the chiral transition in the presence 
of an axial anomaly was argued @ to be second-order and 
in the universality class of the 0(4)-spin model. For finite 
quark masses, due to the explicit breaking of the chiral 
symmetry, the second-order transition is most likely re- 
placed by a rapid crossover. Moreover, arguments based 
on effective models @, i, |, 0, 1, i M, [Li M GJ, M 
indicate that at large n and moderate temperatures the 
transition along the phase boundary is first-order. 

The different nature of the transition at high and low 
/i suggests that the QCD phase diagram exhibits at least 
one critical endpoint (CEP), where the first-order phase 
transition line ends in a second-order critical point, fol- 
lowed by a crossover region [l4|, SH, [TH, QjJ . It is expected 
that the static critical properties of the second-order chi- 
ral endpoint are governed by the universality class of 
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the 3D Ising model [18| ■ The universal properties of the 
phase diagram have recently been studied extensively in 
an effective theory within the functional renormalization 
group (FRG) approach, e.g. @, [H, M, H2 . 

Although there is still no direct proof for the exis- 
tence of the CEP in the QCD phase diagram, many phe- 
nomenological models predict such a point. However, its 
location is still strongly model dependent [3]. Lattice 
QCD simulations also provide indications for a possible 
CEP in the QCD phase diagram Q. 

The critical behavior and the position of the CEP can 
be identified by means of observables sensitive to the sin- 
gular part of the free energy. Of particular interest in this 
context are observables which reflect the fluctuations of 
conserved charges 

H ! i S M HI HI- For systems 
in thermodynamic equilibrium, both baryon number and 
charge density fluctuations diverge at the CEP [lj], in- 
dependently of the value of the quark mass. However, a 
similar divergence is expected also at a first-order phase 
transition in non-equilibrium situations, where the spin- 
odal instabilities are reached [24| . 

The hydrodynamic expansion of an ideal fluid follows 
trajectories of constant entropy, the so called isentropes. 
Due to baryon-number conservation, such trajectories 
correspond to contours of constant entropy per baryon 
s/n in the temperature-chemical potential plane. In 
Ref. [25[ it was pointed out that an expanding system, 
which follows an isentropic trajectory crossing a first- 
order phase transition between the quark-gluon plasma 
and the hadronic phase, is focused towards the CEP, if 
s/n, at a given point {n,T) on the phase boundary, is 
larger in the quark-gluon plasma. This implies that a 
larger range of initial conditions will end up in the vicin- 
ity of the CEP. 

Recently it was argued that the CEP acts as an at- 
tractor for isentropic trajectories also on the crossover 
side, i.e., for values of the chemical potential n smaller 
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than the value at the endpoint [26|, [2lJ- This feature 
could potentially be used to experimentally verify the 
existence of the CEP in ultrarelativistic nucleus- nucleus 
collisions [H, [2§| . 

We note that the relevance of the iscntropic trajecto- 
ries for relativistic nucleus- nucleus collisions rest on the 
assumption that ideal hydrodynamics provides a good 
approximation to the true expansion dynamics and hence 
that dissipative effects can be neglected. Indeed, recent 
analyses of relativistic heavy-ion collider (RHIC) data, 
show that the QGP, created in nucleus- nucleus collisions, 
behaves as an almost ideal fluid with a very small T)/s, the 
ratio of shear viscosity to entropy density (30l | This im- 
plies, that dissipative processes in the QGP are strongly 
suppressed, leading to an essentially isentropic expansion 
towards the phase boundary [3lj . 

However, this is not the case at the CEP [H] where 
dynamical scaling implies that both, the shear and bulk 
viscosities, diverge. # : Hence, close to the CEP the ex- 
pansion is most likely not isentropic and consequently 
the relevance of isentropic trajectories is questionable. 
Closely related to this issue is the critical slowing down 
of long-wavelength fluctuations close to a second-order 
phase transition, which implies that the equilibration 
time of such a system diverges at the critical endpoint. 
It follows that as the CEP is approached, it becomes in- 
creasingly unlikely that an expanding system remains in 
equilibrium, i.e., expands isentropically. However, let us 
for the moment ignore this problem, and assume that the 
expansion of the system is sufficiently slow, so that local 
thermodynamic equilibrium is maintained everywhere. 

Under this assumption, the trajectory of the fireball 
produced in a nucleus-nucleus collision follows an isen- 
tropic trajectory. The corresponding value of s/n is a 
function of the collision energy. Thus, in such a sys- 
tem, the isentropes in the (T, /i)-plane encode important 
information on the expansion dynamics. Consequently, 
possible modifications of the isentropes near the CEP 
influences the evolution of the fireball in a characteris- 
tic way, which could provide an experimental signature 
for the existence of the endpoint. The search for the 
CEP is one of the objectives of planned nucleus- nucleus 
collision experiments at RHIC/BNL J3(| and of future 
experiments at FAIR/GSI [37[. 

In this work we investigate the properties of isentropic 
trajectories near the CEP within an effective chiral the- 
ory. We examine the model dependence of the focus- 
ing on the crossover side of the CEP [2(| both within 
the mean-field approach and in a functional renormal- 
ization group (FRG) analysis, where fluctuations and 
nonperturbative effects are properly accounted for. The 
FRG method has been successfully employed to describe 
a broad range of critical phenomena [f| [01 [2(| H3, [H[ . 



Consequently, this method is ideally suited for exploring 
the effect of fluctuations on the isentropic trajectories. 
Note that, due to the lack of gluonic degrees of freedom 
in the model, the entropy density is much below that of 
QCD matter at the same temperature and chemical po- 
tential. Nevertheless, since the model is most likely in the 
universality class of two flavor QCD, it can provide guid- 
ance on universal properties like e.g. critical exponents 
and the conjectured focusing of the isentropes towards 
the CEP. 

The paper is organized as follows. In the next sec- 
tion we introduce the effective quark-meson model and 
the FRG approach. In Sec. [3] we discuss the mean-field 
thermodynamics, while in Sec. 2] we present the proper- 
ties of the isentropic trajectories near the chiral phase 
transition. Finally, we summarize our results in Sec. [5] 



2. FLOW EQUATION FOR A CHIRAL 
EFFECTIVE THEORY 

The quark-meson model is a low-energy effective the- 
ory, which incorporates the chiral symmetry of QCD. For 
two quark flavors and SU (3) c color symmetry the model 
Lagrangian reads 

1 2 

C = tp [ip - g(cr + ij 5 f ■ tt)} tp + - (d^) - U (a, if), 

(2.1) 

where <fi = (a, if) is the 0(4)-representation for the 
isoscalar a- and the isovector 7?-mesons. The two-flavor 
quark field tp couples to mesons via the flavor-blind 
Yukawa coupling constant g. The purely mesonic po- 
tential is given by 

U(a 7 n) = ^m 2 ^ + ^(^) 2 -ha , (2.2) 

where A > is the mesonic self-coupling. For a negative 
to 2 , the vacuum exhibits spontaneous chiral symmetry 
breaking, where the chiral SU(2)l x SU(2)r symmetry 
of the Lagrangian is broken down to the SU (2) vector 
symmetry. At high temperatures and/or densities the 
spontaneously broken chiral symmetry is restored. The 
external field h, which is related to the current quark 
masses, breaks the chiral symmetry explicitly. 

In the chiral quark-meson model, explicit gluonic de- 
grees of freedom are missing. The effect of gluons is to 
a certain extent implicitly included in the coupling con- 
stants. However, recent extensions of the quark-meson 
model in this context (e.g. by including the Polyakov 
loop) are possible and currently under investigation. 



2.1. FRG equation for the effective potential 



* x An increase of the bulk viscosity near the CEP [33| is suggested 
by model calculations I'{4| and LGT studies |35|| . 



In this section we present the flow equation for the 
thermodynamic potential Q(T,fi) in the FRG approach. 
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Within this scheme, we then compute the thermody- 
namic potential as well as the entropy and baryon num- 
ber densities near the chiral phase transition, including 
the effects of fluctuations. The FRG method yields the 
average effective action Tk at the momentum scale k. 
This scale is introduced in the flow equation via regula- 
tor functions, which act as mass terms in the propagators. 
The effect of the regulators is to suppress the propagation 
of particles with momenta smaller than k. In the infrared 
limit (k = 0) fluctuations at all wavelengths have been 
integrated out and Tk=o is the full effective action. 

The FRG flow equation smoothly interpolates the 
physics between the ultraviolet (UV) A and the IR scale 
[39l . [iol . l4ll ] . For the quark-meson model it contains two 
terms (t = ln(fc/A)) 



d t T k [<f>,*] = ^Tr 



-Ti- 



ft R B ,k (?% [*,*] + Rbx) 



dtRF,k{T ( p[[<f,^} + R F , l 



.(2.3) 



The first term in Eq. (|2.3j) represents the bosonic flow 
with the regulator R B ,k while the second part stands for 
the fermionic contribution with the regulator Rf.u- The 
bosonic and fermionic fields are denoted by $ and ^Z, 
respectively. The full bosonic (fermionic) inverse prop- 



agator includes the term T 



(2) 

B(F).k' 



derivative of Tk with respect to the corresponding fields 
$ or 

As initial condition for the flow equation (|2.3p one 
chooses a bare effective action at the UV scale A. Con- 
sequently, at the momentum scale A, coincides with 
the classical action S, i.e., Tk=\ = S. During the evolu- 
tion, the initial bare action is rcnormalizcd and finally, at 
k = 0, corresponds to the full effective action Tk=o = T- 



2.2. Leading order derivative expansion 



The functional flow equation (|2.3j) is exact and is equiv- 
alent to an infinite tower of coupled partial differential 
equations for n-point functions (n > 2). In order to 
solve this equation, a suitable approximation scheme is 
required. 

In the leading order (LO) derivative expansion the 
FRG flow equation reduces to an ordinary non-linear dif- 
ferential equation [H, [43| . This can be solved either di- 
rectly, for the effective potential on a momentum-space 
grid or by expanding it in powers of the order parameter 
around a minimum. In this letter we use both methods 
and compare the results. For the quark-meson model, in 
four Euclidean space-time dimensions, the LO derivative 
expansion yields [H, IH| 



r fc [^,V>,</>] = / d 4 x [-0 (p + g (a + if ■ n^))^ 



+ \{d^f + U k (p) 



(2.4) 



where Uk(p) is the scale-dependent effective potential, 
cf. Eq. (|2.ip . In this expression we have introduced 
a symmetric field variable p defined by p ~ \<$ 2 = 
\ {a 2 + 7? 2 ) . For a uniform system it is convenient to 
deal with the thermodynamic potential density Q k = 
(T/V)T k = Uk{p = Po,k), where the scale-dependent 
minimum of the potential is labelled by po.fc- In this 
exploratory work, we neglect the scale evolution of the 
Yukawa coupling g. This approximation does not signif- 
icantly affect the critical properties [l^, HH, [H| . 

In thermal equilibrium the extension of the flow equa- 
tion to finite T and p is done within the Matsubara for- 
malism: the integration over pq is replaced by a sum- 
mation over the corresponding Matsubara frequencies: 
Po — > 2mtT for bosons and po — ► (2n + 1)ttT for fermions. 
A finite chemical potential is introduced in the fermionic 
part of the Lagrangian by shifting the derivative with re- 
spect to the Euclidean time r, d T — > d T — p. We assume 
S'C//(2)-symmctry and set p = p u = pd- 

We employ the optimized regulators [HI, [46| for 
fermions 



RfAv) = 



Wo 




the second functional and bosons 



fl B ,fe(q) = (k 2 — q 2 ) (k 2 - q 2 



lb(fc 2 -q 2 )(2.5) 



(2.6) 



where qo = qo + ip- With these, the momentum integra- 
tion and Matsubara summation in the RG flow equation 
can be done analytically, resulting in the following flow 
equation for the thermodynamic potential density Q k 



<9 fe f2fc (T,p;p ,k) 



127T 5 



l + 2n B (E 7r ) 



l + 2n B (E a ) _ ^ 1 - n+(E q ) - n F (E q ) 

En E n 



.(2.7) 



In Eq. (HU) n B {x) = [e x ' T - l] 
Einstein and rip(x) = [ 



= |p(zTM)/T 



M 2 

TT,a,q 



denotes the Bose- 

l] 1 the Fermi-Dirac 
distribution functions for bosons, quarks and antiquarks, 
respectively. The single-particle energies of pions, sig- 

mas and quarks are given by E^^^ = ^ti 2 

where the masses are scale dependent: M 2 = Cl' k , M 2 = 
Cl' k + 2po,fc^fe and M 2 = 2g 2 p ,fc- The prime on the po- 
tential denotes the derivative of Clk = ilk + h \[2~p with 
respect to p evaluated at the minimum po,k- The quark 
degeneracy factor is v q = 2NfN c = 12. 

Although the RG flow equation (|2.7|) looks rather in- 
nocuous, it is quite powerful due to non-linearity implied 
by the self-consistent determination of the single-particle 
energies. This scheme very efficiently accounts for long- 
range fluctuations and nonperturbative dynamics near 
the chiral phase transition (4L|, 5H . 



4 



2.3. Solving the FRG flow equation 



In order to solve the flow equation (|2.7[) , we employ 
two distinct methods. The first one is the grid method, 
where the potential is discrctized on a one-dimensional 
Pi grid. This leads to a set of coupled flow equations for 
the scale-dependent potential f2fc(T, p; pi) at each grid 
point pi. Using the bare potential as initial condition at 
a given UV cutoff A, the scale evolution of the thermo- 
dynamic potential is obtained by finding the minimum 
with respect to variations of the field p at each scale. For 
details concerning the numerical implementation we refer 
toRef. (13. 

The second method is based on a Taylor-expansion, 
up to a maximum power N, of £l k around the scale- 
dependent running minimum p k = °o fc/^ 



The flow equations (|2.11|) - (|2.14|) arc solved numerically 
starting from the initial conditions for the coefficients 
a i,A (i = 0,1,2,3) at the ultraviolet scale A. The model 
parameters of the Lagrangian (|2.1| are chosen such that 
at the scale k = 0, the chiral symmetry is spontaneously 
broken in vacuum with a pion mass of M n = 138 MeV 
and a finite expectation value of the sigma field (a) = <jq , 
which is identified with the pion decay constant f v = 93 
MeV. This, together with our choice for the sigma mass 
in vacuum M a = 670 MeV and the constituent quark 
mass M q = 335 MeV, are obtained with h = 1.771 x 10 6 
MeV 3 , g = 3.6, A = 950 MeV, a hA = (582MeV) 2 , a 2 , A = 
35.2 and = 0. The corresponding values for the 
parameters of the Lagrangian (|2.1|) are A = o.2,a/2 and 
m 2 = «i,a — Actq a , where cto,a = h/ai^ ~ 5.2 MeV is 
the starting value for the scalar condensate. 



N 

&k(T,p;p)= V -H±-(p- po k y 



(2. 



rn—0 



The coefficients a m ^k are functions of the scale k and of 
the temperature and chemical potential. 

On the one hand, the expansion method yields the po- 
tential only in a limited range around the minimum, due 
to the finite convergence radius of the Taylor expansion. 
It is therefore difficult to describe a first-order phase tran- 
sition, where two local potential minima have to be con- 
sidered. On the other hand, the advantage of the Taylor 
expansion is its simplicity; only A^+l coupled differential 
equations have to be solved. 

The minimum of the thermodynamic potential is de- 
termined by the stationarity condition 



da 



da 



-h = 0. 



(2.9) 



In the Taylor expansion scheme this yields a relation be- 
tween the coupling a\ and the expectation value of the 
scalar field ao 



ai,k = h/ao,k ■ 



(2.10) 



We truncate the expansion up to the third order, TV = 3, 
and obtain the following set of flow equations: 



dk.ao. 



dkPo,k 



y/2po,k 



d kPo ,k + d k n k , (2.11) 



(o 12") 

h/(2p ^ 2 + a 2 , k ' 1 ' ' 

dka 2 ,k = 03, fc d k po.k + d k , (2.13) 

d k a 3 , k - d k O,'^ . (2.14) 



The meson masses are then given by 



Mik = 



Mi k = 



y/2>P0,k 



+ 2p 0:k a 2 ,k, (2.15) 



2.4. Isentropic thermodynamics 



The solution of the flow equation (|2.7[) yields the pres- 
sure p(T,p) = — flk{T, p\ Po,k)\ k=0 asafunction ofTand 
p. The relevant observables, the entropy density s{T,p) 
and the quark-number density n(T, p) are obtained from 
the pressure as derivatives with respect to temperature 
and quark chemical potential 



dp(T, p) _ d [a ,fe - ha 0tk ] 



dT dT 
dp(T, p) d [a ,fc - ha 0tk ] 



, (2.16) 



dp 



dp 



k=0 



k=0 



(2.17) 



The isentropic trajectories in the (T, ^i)-plane are then 
obtained as contours of constant entropy per baryon s/n. 



3. MEAN-FIELD APPROXIMATION 

In order to examine the influence of fluctuations on the 
thermodynamics, in particular on the isentropic trajec- 
tories, we compare the FRG results with those obtained 
in the mean-field (MF) approximation, where quantum 
and thermal fluctuations are neglected. 

The partition function of the quark-meson model 
can be formulated as a path integral over meson and 
quark/antiquark fields in Euclidean space-time. In the 
MF approximation, the meson fields in the action are re- 
placed by their expectation values. The integration over 
the quark/antiquark fields yields the fermionic determi- 
nant. The resulting thermodynamic potential of the chi- 
ral quark-meson model is of the form 0, [48| 

n(T,p; (a), (7?)) = n Qg (T,p; (a)) + U((a), (7?)) (3.1) 



and the quark mass by M 2 k = 2 g 2 po.fe- 



with the meson potential U of Eq. (|2.2p . The quark/anti- 
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quark contribution is given by 



d 3 k 



{\n(l -n F (E q )) 



2tt) 3 

ln(l-n+(E,))-^} ,(3.2) 



where the last term is the divergent vacuum contribution. 
As will be shown below, this term influences the shape 
of the iscntropic trajectories near the phase boundary. 




50 100 150 200 250 300 350 
H [MeV] 



FIG. 1: Isentropes computed in the mean-field approxima- 
tion to the quark-meson model, neglecting the vacuum term 
in Eq. (|3.2p . The values of s/n are indicated at each isen- 
tropic trajectory. The chiral phase boundary, composed of a 
crossover and a first-order transition, is indicated by a broken 
and a full line, respectively. The bullet on the phase boundary 
indicates the position of the CEP in the MF approximation 
to this model. 

In the MF approximation the expectation value (a) 
is determined by the corresponding classical equation of 
motion, the gap equation. This is obtained by minimizing 
the thermodynamic potential in the er-direction 



on 



= 



(3.3) 



The solution of the gap equation determines the T- and 
/i dependence of the chiral order parameter (a)(T, /j) and 
the constituent quark mass M q = g{a). The expectation 
values of the pion fields (7?) vanish. We have chosen the 
model parameters in the MF analysis so as to reproduce 
the same vacuum physics as in the FRG approach. 

In Fig. Q] and [2] we show the isentropes obtained in 
the mean-field approximation with and without the vac- 
uum contribution. The divergence of the vacuum term 
is regularized by a momentum cutoff A = 583 MeV. The 
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FIG. 2: The same as Fig. [T] but including the vacuum term 
in Eq. (pHJ) . 



parameters of the model were in both cases chosen so as 
to reproduce the vacuum physics discussed in sectior l2~3l 
We note that the position of the CEP is shifted to smaller 
temperature and larger values of the chemical potential, 
namely from (T c , fi c ) = (130, 147) to (51, 316) MeV, when 
the vacuum contribution is included. Furthermore, the 
crossover temperature at fi = is increased from ap- 
proximately 150 to almost 210 MeV, while the critical 
chemical potential at T = is changed by less than 10 
%. 

Our results also show that the kink in the isentropic 
trajectories found at the phase boundary in the calcula- 
tion, where the vacuum contribution is neglected, is to a 
large extent an artifact of the approximation. The ther- 
mal part of (|3.2[) yields a term of the form logM?. 
With such a term in the effective potential, the phase 
transition in the chiral limit is first order for all densi- 
ties, from fi = and T w 200 MeV to T = and fi w 300 
MeV as recently found in Ref. [f|. Due to the non- 
zero latent heat, the isentropic trajectories are discontin- 
uous at the phase boundary in the T — fj, plane. When 
a sufficiently strong explicit symmetry breaking term is 
introduced, the first order transition is smoothened into 
a crossover and the discontinuity in the isentropes shows 
up as kinks at the phase boundary. However, when the 
vacuum term is included, the logarithmic term is can- 
celled and the transition in the chiral limit is second or- 
der in flavor SU(2). Consequently, in the chiral limit the 
isentropes are continuous and when the explicit symme- 
try breaking term is turned on the isentropic trajectories 
remain smooth, as seen in Fig. [2] 

These considerations offer a plausible explanation for 
the different behavior of the isentropes near the phase 
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boundary in the flavor SU(2) versions of the quark- 
meson and Nambu-Jona-Lasinio (NJL)) models found in 
Ref. [HI]. In the quark-meson model the vacuum term 
was dropped, while in the NJL model it must be included, 
since the model otherwise does not yield a state with 
spontaneously broken chiral symmetry. Consequently, 
one expects a kink at the transition in the quark-meson 
model and smooth isentropes in the NJL model, in agree- 
ment with the results of Ref. [48| . We conclude that the 
strong kink structure in the crossover region, which is 
obtained by neglecting the quark vacuum loop, is un- 
physical. 

In two recent papers, the isentropes have been com- 
puted in effective chiral models including the coupling 
to the Polyakov loop in the mean-field approxima- 
tion [Hol . \EM- Also in this case, the singular term 
cancels between the thermal and vacuum contributions. 
Hence, one expects a kink in the quark-meson model with 
Polyakov loop and smooth trajectories in the Polyakov- 
NJL (PNJL) model. This is confirmed on a qualitative 
level by the calculation of Ref. [5(| and similar trajecto- 
ries have been obtained in the flavor SU(3) version of the 
PNJL model in Ref. [5l[ . The cause of a less pronounced 
structure found in the isentropes obtained in the PNJL 
model as the CEP is approached, remains unclear. A 
possible connection with the first order transition in the 
chiral limit of the flavor SU (3) model [52| is uncertain, 
since similar trajectories are obtained in the two flavor 
PNJL model, where in the same limit the transition is 
second order [Hoj . 

At the first-order transition to the broken symmetry 
phase, the isentropes are deflected to larger \x and smaller 
T, away from the CEP, because the entropy per baryon in 
the co-existence region is, in chiral effective models lack- 
ing gluonic degrees of freedom, smaller in the symmetric 
phase. This may be different if the chiral and deconfinc- 
ment transitions coincide at finite n, since in QCD the 
gluons, are expected to contribute a major part of the 
entropy per baryon. However, if the so called quarkyonic 
phase [53J is realized, the isentropes at the first-order 
transition from hadronic to quarkyonic matter may be 
similar to those found in the quark-meson model. 



4. ISENTROPIC TRAJECTORIES IN THE RG 
APPROACH 

As indicated in the introduction, the influence of fluc- 
tuations on the isentropic trajectories at the chiral tran- 
sition and in particular at the critical endpoint of QCD is 
of high current interest [H, [H, [2^1 . We address this prob- 
lem by comparing the isentropic trajectories obtained us- 
ing the FRG approach with those obtained in the MF 
approximation in the preceding section. 

Using the thermodynamic potentials obtained in the 
FRG approach, we determine the phase diagram of the 
model [54j. The crossover line corresponds to a maxi- 
mum of the chiral susceptibility. In Fig. [3] we show the 



resulting phase diagram together with the corresponding 
isentropes. As in the MF approach, the phase diagram 
has a generic structure and exhibits a critical endpoint, 
which separates the first-order phase transition from the 
rapid crossover. At the CEP, the transition is of second- 
order and belongs to the universality class of the Ising 
model in three dimensions. 

A comparison with the MF results shows that the FRG 
phase diagram is very similar to that obtained in the 
mean-field calculation with the vacuum term included. 
In particular, the position of the critical endpoint, the 
crossover temperature at fi = and the critical chem- 
ical potential at T = are almost identical. Thus, 
for the quark-meson model, the effect of fluctuations on 
the phase diagram are fairly small, once the fermion 
vacuum term is taken into account. With our choice 
of model parameters the FRG approach yields approx- 
imately (T c , (j, c ) = (14, 328) MeV within the grid method 
(FRGg r id), while using the Taylor expansion method 
(FRGxayior) we find (T c , (i c ) = (64, 336) MeV, as shown 
in Fig. [31 Thus, the location of the CEP is not only 
model dependent [l4|, but also depends on the way the 
dynamics is implemented in the renormalization group 
flow equation. The two approaches differ in the treat- 
ment of higher order, irrelevant operators. This does not 
affect universal quantities, like critical exponents, but it 
can affect non-universal quantities, like the location of 
the critical endpoint. 

A shift of the CEP to larger values of the chemical po- 
tential and smaller temperatures was also found in the 
proper-time renormalization group approach relative to 
the mean- field solution in Ref. [5|. In view of our re- 
sults, it is likely that a major part of this shift is due to 
the fermion vacuum term, which was not included in the 
mean-field approximation. 

We also show the isentropic trajectories in the (T, /x)- 
plane, obtained in the FRGTayior and in the FRG gr id ap- 
proaches. The isentropes show a smooth behavior every- 
where, and agree qualitatively with the mean-field tra- 
jectories shown in Fig. [2l i.e., those obtained with the 
vacuum term. It is reassuring that the FRGxayior and 
the FRGgrid methods lead to a very similar isentropic 
trajectories. Some deviations at large /x, seen in Fig. [3J 
indicate limitations of the Taylor expansion method due 
to the truncation of the series (|2.8[) at third order. Re- 
cent LGT results [55| on the isentropic trajectories near 
the QCD crossover line show qualitatively a very similar 
smooth behavior as obtained in our model with the FRG 
approach and shown in Fig. [3J and [2] 

A comparison of the FRG and MF results shows no 
qualitative change of the isentropic trajectories in the 
vicinity of the CEP. In particular, the isentropes remain 
smooth also when the effect of long-wavelength fluctu- 
ations is consistently included and show no sign of fo- 
cusing towards the CEP. This is clearly illustrated in 
Fig. [H which shows a blow-up of the FRG isentropes 
in the vicinity of the CEP. Since the quark-meson model 
is in the same universality class as QCD, our results do 
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yu [MeV] 

FIG. 3: Isentropes calculated in the quark-meson model us- 
ing the FRGTayior and FRG gr id methods (see text). The isen- 
tropes computed within the FRGTayior approach are shown 
as solid lines whereas those obtained within the grid method 
FRGgrid are indicated by dashed lines. The s/n-ratio of each 
isentrope is indicated at each contour. The phase boundary, 
with the CEP, obtained in the FRG approach, is indicated as 
in Fig. [1] The CEP shown is that obtained using the Taylor 
expansion method FRGTayior • 



not support the universality of the focusing phenomenon 
conjectured in Refs. [26l.l28j]. We stress that the RG treat- 
ment of fluctuations employed here, reproduces the Z(2) 
universal scaling of the relevant physical observables at 
the CEP 0. 

The fact that the focusing effect is not universal can be 
understood in general terms. The point is that the en- 
tropy and the baryon density are both obtained as first 
derivatives of the thermodynamic potential Q, which re- 
main finite at the CEP, since only second- and higher- 
order derivatives diverge. Consequently, the singular 
part of the entropy per baryon does not diverge at the 
CEP and hence is not guaranteed to dominate over the 
regular background contribution. It follows that the isen- 
tropic trajectories are not universal, since they depend on 
the relative strength of the universal singular part and 
the non-universal background. In other words, the char- 
acteristic shape of the isentropes in the vicinity of the 
CEP can vary from model to model, even though they 
belong to the same universality class. The model con- 
structed in Refs. 0, [28| yields focusing of the isentropes 
towards the CEP because the singular part of the thermo- 
dynamic potential is chosen by hand to be very large. In 
chiral models, where the critical region around the CEP 
and around the 0(4) transition line is small 0, HH, an( i 
consequently the relative strength of the singular part of 
f2 is small, it is unlikely that the focusing effect of the 



isentropic trajectories reported in (26j can be observed. 




320 322 324 326 328 330 
H [MeV] 

FIG. 4: Isentropes as in Fig.[3]but in a narrower region around 
the CEP from the side of the crossover transition (dashed 
line) . The CEP is indicated as a black point at the right edge 
of the figure. 

For completeness we note that the fact that this effect 
is not universal does not exclude the possibility that in 
QCD the relative strength of the singular part is large. 
If this is the case, focusing may then potentially be rel- 
evant for nucleus-nucleus collision experiments (28l . |29| . 
However, then a detailed study of the equilibration of 
long-range fluctuations in an expanding system, along the 
lines of Ref. [I?) is needed in order to decide whether the 
isentropic trajectories are relevant or not, as discussed in 
the introduction. 



5. SUMMARY 

The isentropic trajectories (contours of fixed entropy 
per baryon) in the QCD phase diagram describe possi- 
ble paths of the hydrodynamic evolution of a thermal 
medium created in nucleus- nucleus collisions. We inves- 
tigated the behavior of the isentropic trajectories within 
the chiral quark-meson model for two-quark flavors. The 
thermodynamics was formulated using functional renor- 
malization group (FRG) techniques and the results were 
compared with two variants of the mean-field (MF) ap- 
proximation, one neglecting and the other one including 
the fermion vacuum term. 

Our studies of the isentropic trajectories near the chiral 
phase transition were motivated by recent findings that 
the chiral critical endpoint (CEP) acts as an attractor 
for the isentropes, leading to a focusing towards the CEP 
[26[. It was argued that the focusing effect would have 
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important phenomenolqgical consequences for nucleus- 
nucleus collisions [H, [291 ]. 

A comparison of the MF and the FRG results for the 
isentropes show that the kink structure in the transi- 
tion region, observed in some MF calculations, is washed 
out when the fermion vacuum contribution is properly 
included. Furthermore, in spite of the fact that the en- 
tropy density and the baryon-number density rise rapidly 
near the crossover transition and also as the CEP is ap- 
proached, the isentropes around the CEP remain very 
smooth. These results raise doubts concerning the focus- 
ing of isentropes and its phenomenological consequences. 

Although, as already mentioned, the isentropic trajec- 
tories of our model differ from those of QCD, the singular 
part of the thermodynamic potential near the CEP is uni- 
versal. Our results show that in the quark-meson model 
the singular contribution to the entropy per baryon near 
the CEP is subdominant and hence that the conjectured 



focusing of the isentropic trajectories towards the CEP 
is not universal. However, the possibility still remains 
that in some systems of the same universality class fo- 
cusing may appear, if accidentally close to the CEP the 
singular term dominates. Thus, we conclude that the ex- 
istence of such an effect in QCD matter is unlikely but 
not completely excluded. 
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